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I. INTRODUCTION 

Quantum phase transitions (QPTs) [1] in spin systems, 
e.g., the XY model % the Lipkin-Meshkov-Glick (LMG) 
model Q, and the Dicke model [4], have aroused much 
interest in recent years. Most of these efforts have ad- 
dressed possible connections of quantum entanglement 
measures, such as the concurrence, the entanglement en- 
tropy, and the negativity ^ith the QPTs in the systems. 
The scaling behavior [2, Q of the entanglement demon- 
strate well the quantum criticality of these systems. On 
the other hand, being related to quantum measurement 
theory and quantum decoherence problems, theoretical 
studies of the bath influence on the dynamic property 
of a central system have also attracted much attention. 
Ref. claimed to find "an intrinsic limit to quantum 
coherence due to spontaneous symmetry breaking" in 
Lieb-Mattis model; the relationship between entangle- 
ment dynamics and paramagnet-ferromagnet phase tran- 
sition was explored in Ref. Q; it was shown in Ref. [8] 
that the Loschmidt echo decay enhanced at the criti- 
cal point can be a signature of QPT in the transverse 
Ising model. In fact, further deeper studies on the the 
relevant issues of open quantum systems not only pro- 
vides us a better understanding of the quantum-classical 
crossover, but also promises important potential applica- 
tions in quantum information processing [9]. 

In this paper, integrating coherently the above two 
interesting topics: quantum phase transition and quan- 
tum open system , we elaborate how the QPT of the 
"bath" influences the dynamics of a central qubit cou- 
pled to it. It is shown that when the bath is in different 
phases, the purity of the central qubit exhibit distinctly 
different behaviors in two different phases. Moreover, it is 
also illustrated that under certain conditions our model is 
equivalent to an anti Jaynes-Cummings (anti J-C) model 

The paper is organized as follows. In Sec. II, we intro- 
duce the LMG model and summarize the main properties 
of this model. In Sec. Ill, we study the dynamic evolu- 
tion of a central qubit coupled to a bath described by 



an isotropic LMG model, which is exactly solvable. In 
Sec. IV, We evaluate the purity of the central qubit in 
the symmetry broken phase and the symmetric phase, re- 
spectively. The QPT of the bath is well indicated by the 
behavior of the purity. In Sec. V, a connection between 
the current model and an anti J-C model is established. 
Section VI presents our summary and conclusion. 



II. LIPKIN-MESHKOV-GLICK (LMG) MODEL 
FOR QUANTUM PHASE TRANSITION 

We consider a central qubit (two-level system) that 
couples to a multi-qubit bath, which is described by the 
LMG model H 
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where erf, a = x,y, z [i = 1,2, ■ ■ ■ N) are the Pauli matri- 
ces of the z-th atom, \/N denotes the coupling strength, 
which is inversely proportional to the atom number N . 
This Hamiltonian contains long-range interactions, i.e., 
every spin in the bath interacts with all the others. In 
the isotropic case, 7 = 1, the Hamiltonian is diagonal in 
the Dicke representation 
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and the ground state of Hb lies in the subspace spanned 
by the Dicke states {\N/2, M),M = -N/2, ■ ■ ■ N/2} 
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The eigenenergy corresponding to \N/2, M) is 2AM^/iV— 
2M — XN/2. Hence, the ground state |G) is A-dependent 
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where 



|G) 



f,f) , (0<A<1), 
"•,/(A)), (A>1) 
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where /(A) is the integer nearest to N/2X. The level 
crossing at A = 1 leads to the occurrence of a QPT. This 
point is also a symmetry breaking point p^ . [3] : when 
< A < 1, the ground state of the bath is unique and 
fully polarized in the magnetic field direction, and thus 
the bath is in a symmetry broken phase; when A > 1, the 
ground state is infinitely degenerate and thus the bath is 
in a symmetric phase. We below elaborate how the dy- 
namic evolution of the purity (a measure of quantum 
coherence) depends on the coupling strength between the 
central spin and the bath; in particular, we observe that 
the purity shows distinctly different behaviors in the two 
phases, which may be used to reveal the QPT point in 
the bath. 



/3 = [N^ - 4(Af + 1)2] - 2(Af + 1) + 1, 



(8) 



C = -AV^ (A^ + 2) - 4M{M + 1). 

A straightforward calculation determines the two eigen- 
values xi and X2 of Hm as 



XI = i[(a + /3) + V(«-/3)2+4C2], 
X2 = i[(a + /3)- V(« -/?)'+ 4C2], 



and the eigenstate j^E'i) corresponding to xi is 
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III. DYNAMICS OF A CENTRAL QUBIT 
COUPLED TO AN ISOTROPIC LMG MODEL 

A spin-bath model is described by the total Hamilto- 
nian H ^ Hb + Hs + Hsb where Hs = -2s^ is 

the free Hamiltonian of the central qubit S. Hsb denotes 
the coupling between S and the bath B. Specifically, the 
total Hamiltonian can be written as 
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where cr", a = x,y,z, are the Pauli operators of the 
central qubit; A' is the coupling strength between the 
central qubit and the bath. In the Dicke representation, 
the above Hamiltonian can be rewritten as 



H 



A 

'n 



[JnI 



N 



Jn^n 



N\ - 2J, 



N 



(6) 



2\'{s+J^ + s^J+)-{2s,), 



where = J^±iJff and s± = s^±is^ are the ladder op- 
erators of the iV-qubit bath and the central qubit, respec- 
tively. For simplicity, we denote the two eigenstate of the 
central qubit as |t)' = |l/2,f/2) and |i) = |l/2,-l/2), 
and s^lt) = IT) /2, s^H) = -\i)/2. In an invariant 
subspace Ti-M of H spanned by the ordered basis vector 
{|7V/2, M) ® It) , |A^/2, li)}, the total Hamilto- 

nian can be expressed as a quasidiagonal matrix with the 
diagonal blocks 



Hm — 



a, C 

C, P 
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where 



c 



xi — a 
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We would like to point out that, in the symmetric phase, 
all xi, X2, a, and b are functions of /(A), i.e., F = F[I (A)] 
with F = xi, X2, a, and b. The dynamic evolution op- 
erator U(t) — eKp[—iHt] in the subspace Ti.M can be 
expressed in terms of xi, X2, a, and b 
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ab{e 
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We wish to mention that the above dynamic evolution 
(fTT|) is valid only for the cases when -A^/2 < M < N/2, 
because {\N/2, M) (g) |t) , |iV/2, M + I) ® |i)} is a two- 
dimensional invariant subspace for these cases. But for 
the case M A^/2, {|7V/2, M)(g) |t)} is a one-dimensional 
invariant subspace, i.e., \N/2, N/2) (E) |t) is an eigen- 
state of the total Hamiltonian ([5]), and its corresponding 
eigenenergy is —{N +1). Thus the dynamic evolution of 
this state \N/2, N/2) (g) |t) is different from Eq. We 
will discuss this point in the next section. 



IV. PURITY OF THE CENTRAL QUBIT AS A 
WITNESS OF QUANTUM PHASE TRANSITION 

Based on the above results, we now solve the 
Schrodinger equation that describes the dynamics of the 
purity of the central qubit. To highlight the influence of 
the QPT of the bath on the coupled central qubit, it is as- 
sumed that the bath and the central qubit are initially in 
the ground state \G) dH) and a pure superposition state 
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C| It) + ci 11), respectively. The evolution of the total 
system (the bath plus the central qubit) is 



|*A.+i(t))=e-'«*|G)®(cT|T)+cx|i)), (12) 
and the reduced density matrix of the central qubit is 



p^(t) =TrB|^'jv+i(t)) {^N+i{t)\, 



(13) 



where Tre means tracing out the degree of freedom of 
the bath. 

The purity P of the central qubit is defined as 



P = TrsiipHm 

= [pfTWP + buW]' + 2|pfx(i)| 



(14) 



which can be used to measure the quantum coherence. 
For a pure state, the purity equals to unity, while for a 
mixed state the purity is less than unity. The decay of 
purity indicates the loss of quantum coherence [lB| ■ 



A. Purity of the central qubit in two phases 

1. Symmetric phase 

When A > 1, the bath is in the symmetric phase and 
/(A) < N/2. As mentioned above, we can apply the 
evolution matrix ([TTjl to obtain the reduced density ma- 
trix {t) of the central qubit with the matrix elements 
defined by 



pfT(i) = ^jM'f{Kt) + \c,\'h{X,tW, 
pfiit) - ^[hfgi\,t) + \c,f^i\,t)]^ 



(15) 



g(A,t)x/(A,t) = |pf^(t)| 



where 

/(A,t) 
5(A,i) 
h(X,t) 
i{X,t) 



+ 2a''6^cos[(a;i - 2:2)*], (16) 
a-y + {h' f + 2{a'f{b'f cos[{x\ - x'^)t\ , 
2(a')'(&')'{l-cos[(a;'i-4)i]}, 
2a^b^{l ~ cos[{xi - X2)t\}, 



and the parameters and b' are defined by 

F' = F[I{\) - 1] with F' = x'l, 4, a', and b' . This 
subtle change from {xi, X2, a,b) to {x'i,X2,a' ,b') is due 
to the fact that \G) (g) |t) and \G) (g) |i) belong to two 
different invariant subspace Hm and Hm-i- Ifj for sim- 
plicity, we assume that the central qubit is initially in 
the superposition state (|t) + \i))/V2, we obtain from 
Eq. (I14p the exact expression of the purity of the central 
qubit 

P ^ hfiX,t) + h{\,t)]' + hziX,t)+giX,t)fil7) 



--g{X,t)f{X,t). 



2. Symmetry broken phase 

When < A < 1, the bath is in the symmetry broken 
phase and the ground state is the fully polarized state 
\G) = \N/2,N/2). \G) (g) It) is an eigenstate of the to- 
tal Hamiltonian ([5]) , and its corresponding eigenenergy is 
— {N + 1). Thus the dynamic evolution of this state is 
exp[i(A^-|- l)t]. While the dynamic evolution of the other 
state \G) (g) II) can be obtained following the way men- 
tioned above Eq.((Tll) with M = N/2 - 1. After a similar 
procedure, if the central qubit is initially prepared in the 
superposition state (|t) + U))/ V^, the exact expression of 
the purity in the symmetry broken phase can be written 



P 



where 



[l + h{X,t)Y + -~g'{X,t) + -g{X,t), (18) 



g{X,t) = a^ + V + 25,%'^cos[{S:i-X2)t], 
h{X,t) = 2a^62{l - cos[(ii - £2)*]}. 

The new parameters xi, X2, a, and b are given by -F = 
F[N/2 - 1] with F = ii, f2, a, and b. 

To obtain the exact result of the purity of the cen- 
tral qubit, we need to know the coupling strength A' [l^ 
between the central qubit and the bath. In different ref- 
erences, this coupling strength is treated differently. For 
example, this coupling strength was assumed to be in- 
verselyproportional to the qubit number of the bath in 
Ref. [6|, while it was assumed to be inversely propor- 
tional to the square root of the spin number in some 
other references 0, [l!]. We hereafter denote the two 
cases with the above two different coupling strengths as 
Cases I and II, respectively. Generally speaking, when 
the central qubit is identical to the qubits of the bath, 
the coupling strength between the central qubit and the 
bath should be equal to the coupling strength between 
the qubits of the bath. This is Case I, and we will elab- 
orate it in next subsection. Besides Case I, Case II will 
also be discussed later. The calculation for Case II is 
the same as that for Case I except that A' in Case I is 
changed to VlVX'. The two different coupling strengths 
will lead to different behaviors of the purity. 



B. The coupling strength inversely proportional to 
the spin number of the bath (Case I) 

Firstly let us consider Case I. Similar to the coupling 
mechanism in Ref. 0|, we assume that the coupling 
strength A' is just the coupling strength X/N between 
any two qubits of the bath. Since the coupling strength 
X/N between qubits in the LMG model is inversely pro- 
portional to the spin number N, the system is extensive. 

Fig. 1 clearly shows that, in the symmetry breaking 
phase, the purity of the central qubit remains as a con- 
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FIG. 1: (Color online) Two different view angles on the dy- 
namic evolution of the purity P as functions of A and t (Case 
I). The QPT (symmetry breaking) at A = 1 is well signatured. 
The purity saturates when becomes large. Here we have 
chosen the qubit number of the environment A*' = 5000. 




symmetry broken phase of the bath can also be verified 
through another approach. The ground state of the LMG 
model in the symmetry broken phase (0 < A < 1) is 
\G) = \N/2,N/2). The direct product of the bath \G) 
and the central qubit c-f |t) + ||) can be expanded in 
the angular momentum coupling representation 



|G)®(ct|T>+C2U)) « 



i(7V + l),i(iV + l))(19) 



+ci \{N-l)^-{N-l) 



where we have used the Clebsch-Gordan (C-G) coefficient 

(20) 



y,y)«IT> 



i(iV + l),i(7V + l) 
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(7V + 1),-(7V-1) 



The total Hamiltonian ^ of the qubit and the bath can 
be rewritten as 



H 



A 



N 



'N+l- 

(21) 

Through the above approximation (|19p . both 
|iV/2, N/2) (g) It) and |iV/2, N/2) (g) |i) are the eigenstate 
of the total Hamiltonian with the eigenenergy —{N + 1) 
and 2X/N — {N — 1) respectively, and then the dynamic 
evolution of the two states are obvious. After a straight- 
forward derivation, the reduced density matrix of the 
qubit is expressed as 



p'it) = icTriT)(Ti + iqru)ai 

-fc,cIe^(^+2)*|T) (il 



(22) 



h.c. 



and the purity of (t) remains as unity by applying Eq. 
(fT4| . It is thus proven that the qubit preserves its quan- 
tum coherence when the bath is in its symmetry broken 
phase (0 < A < 1). 



FIG. 2: (Color online) Dynamic evolution of the purity P as 
a function of time t for different A. In the symmetry broken 
phase (0 < A < 1), P remains as a constant unity. In the 
symmetric phase (A > 1), P ranges from 0.5 to 1, and the 
period of the oscillation decreases as A increases. The curves 
with different colors represent A = 1.01, A = 1.1, A — 1.3, 
A = 2, and A = 5 (AT = 5000). 

stant unity, i.e., the qubit preserves its quantum coher- 
ence all the time. The decay of purity becomes van- 
ishingly small in the thermodynamic limit. While in the 
symmetric phase, the purity varies periodically, as shown 
in Fig. 2. 

The constant purity P = 1 of the central qubit for the 



C. The coupling strength inversely proportional to 
the square root of the spin number of the bath 
(Case II) 

We now turn to Case II 0, ■ Being different from 
Case I, the purity of the central qubit does not preserve 
its coherence when the bath is in the symmetry broken 
phase (see Fig. 3), although the purity varies also period- 
ically in the symmetric phase (Fig. 4). However, both the 
range and the pattern of its time-dependance in the two 
phases are different from those in Case I. We also remark 
that the purity P saturates in the symmetry broken phase 
when N increases, just like that in the symmetric phase 
of Case I; while the dynamic behavior of the purity in the 



symmetric phase depends on N with the period being in- 
versely proportional to \/]V approximately, as analyzed 
later. 
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FIG. 3: (Color online) Dynamic evolution of the purity P 
as functions of A and t (Case II) in the symmetric phase and 
symmetry broken phase. Clearly the purity varies in distinctly 
different manners in the two phases, which may be considered 
as an indication of QPT at the critical point A = 1. The purity 
reaches a steady state in the symmetry broken phase, while 
the pattern of the purity will always change with A*' in the 
symmetric phase. Here we choose A'' = 1000. 



V. EQUIVALENCE TO AN ANTI 
JAYNES-CUMMINGS MODEL 

In this section we show the equivalence between the 
above model and an anti J-C model with an intensity- 
dependent coupling strength 19]. With this observation, 
we can exactly solve the dynamical equation about time 
evolution. 

A. Symmetry broken phase 

When the bath is in the symmetry broken phase, our 
model may be recast into an anti J-C model. Actually 
the equivalence between the LMG model and Dike model 
was just studied recently (20j . 

In the symmetry broken phase (0 < A < 1), the ground 
state \N/2, N/2) of the bath corresponds to a low exci- 
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FIG. 4: (Color online) Dynamic evolution of the purity P as 
a function of time t for different A. In the symmetry breaking 
phase (0 < A < 1), the dynamic behavior of P saturates 
when A'^ increases, which ranges from about 0.88 to 1 with the 
period about 1.4. In the symmetric phase (A > 1), P ranges 
from 0.5 to 1, and the period of the oscillation decreases as A 
increases. The curves of different colors represent A = 1.0001, 
A = 1.0003, A = 1.002, A = 1.02 (TV = 1000). 

tation Fock state |0) after the H-P transformation. The 
mean photon number n — (^d^d^ = 0. Hence we can 
directly expand the Holstein-PrimakofF (H-P) transfor- 
mation f^H to the first-order [22] 

4 = V]vd,j^ = (j+)t 

= N/2-dU, 
and the Hamiltonian ([6]) can be rewritten as 

H = 2{1- X)d^d -N- 2\'y/N{s+d^ + s^d) - 2s,. (23) 

Let us recall that the anti J-C Hamiltonian can be rewrit- 
ten as 

Hajc = vd'^d - k{cr+d'^ + a^d) + ^^^cTz, (24) 

where d^ and d are the creation and annihilation opera- 
tors of the single-mode quantized field with the frequency 
!/; — fc is the coupling strength between the field and the 
two-level system; uj is the level spacing between the two- 
level system; a+ — {ax + iay)/2 and ct_ — {ax — iay)/2. 
Hence, the model described by Eq. (|23p is an anti J-C 
model dm) with = 2(1 - A), w = -2, and k = 2\'^/N. 
We now illustrate that the boson mode characterized by d 
and may be mapped from the collective spin J in a low 
excitation limit. Note that the different mapping ways 
depend on the phases of the bath because the bosoniza- 
tion of collective spin is essentially a mean field approach 
based on choice of the order parameter. 

The solution |^'(t)) of the Schrodinger equation 
ihdt \^{t)) — Hajc \^{i)) can be expressed as 

l*W> = I][cT,n+i W IT) ®\n+l)+ cx,„(i) li) ® \n)], 

n=0 

(25) 
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where n = (^(fid^ is the mean "photon" number. A 
straightforward calculation determines the probability 
amplitudes [llj 



CT,n+l(i) = {ct,„+i(0) 



COS (flnt) — ijr- sin(r2„t) 



(26) 



+i — \/n+ lc|,„(0) sin (ri„i)} exp ( j i 



Ci,n(i) = {ci,n(0) 



A 

COS + i— — sin 



(27) 



+i + lc|,„+i (0) sin ($!„<)} exp ( 



where A = i/ + and ^ y^(A/2)2 + A:2(n + 1). 

After a routine calculation, we obtain the purity of the 
central qubit (with an initial state (|t) + |i))/V2) in the 
symmetry broken phase (0 < A < 1) (see Appendix A) 



P 



1 + cos^ {Q.Qt) + ( ^ ) sin^ {VLot) 



(28) 



Aj sin^ (riot) + 2 A") sin2(r!oi) 



In Case I (A' = A/A^), the coupling strength k — 
2\/\/N is inversely proportional to the square root of 
the qubit number N . In the large N limit, 



N 



I.e., 





/ k ^ 


4fc2 




lim 






= 0, 






/ ^ A2" 




A ^ 






lim 






= 1. 


Af^oo \ 




1 A2 



(29) 



(30) 



Hence from Eq. (l28l) . we have P = 1 in large N limit. 
This analytical analysis agrees well with Eq. (|^^ and 
Fig. 1. 

Generally speaking, the quantum coherence (measured 
by purity) of a quantum open system would be dissipated 
by its bath. But when the coupling strength between the 
system and its bath becomes vanishingly small, the sys- 
tem and bath becomes decoupled, and then the system 
will preserve all its coherence (remains in a pure state or 
the purity remains to be unity) during the dynamic evolu- 
tion. In the thermodynamic limit, the coupling strength 
between the two level atom and the "radiation field" 
becomes vanishingly small, i.e., the radiation field and 
the central qubit are decoupled. Thus, the central qubit 
evolves under the free Hamiltonian Hs = —25^, which 
preserves quantum coherence of the central qubit. 

In Case II (A' — \/^/N), the coupling strength k = 2A. 
Even in the thermodynamic limit, the interaction Hamil- 
tonian does not vanish {2X'^/N = 2A ^ 0). This is why 



the purity of the central qubit does not preserve in Case 
II even when the bath is in the symmetry broken phase. 
In this case, the purity ((28|) can be simplified to 



hm P^- 



— sm^iVbXt) - ^ sin^iv^Xt) + 4 
25 5 



(31) 

which is independent of N, as we observed numerically 
in Sec. IV. 



B. Symmetric phase 

In the above discussion, the system and the bath is 
reduced into the anti J-C model with a normal coupling, 
i.e., the coupling does not depend on the mean "photon" 
number. We now show that the system and the bath can 
be reduced into the anti J-C model with an intensity- 
dependent coupling when the bath is in the symmetric 
phase. In the symmetric phase (A > 1), however, the 
ground state \N/2, 1 (X)) of the bath is no longer a low 
excitation state after the H-P transformation. The mean 
"photon" number 



1 



(32) 



is of the same order of N. By applying H-P transforma- 
tion the Hamiltonian (16]) can be rewritten as 



2A 

H = —{dUf + 2{l-X)dU-N 



(33) 



-2X'{s+dWN - dU + .s_ \/N - dUd) - 2s,_. 



The model described by this Hamiltonian (|33p is an 
intensity-dependent coupling anti J-C model with a Kerr- 
effect term 2X{d^d)'^/N, which can be analytically diag- 
onalized as well. After a similar derivation to that in the 
symmetry broken phase, we obtain the solution |$(t)) of 
the Schrodinger equation ihdt \ ^{t)) = H |$(t)). 



Mt)) - EK,n+i(*) It) (^\n+l)+ c\Jt) \i) (g> In)], 

(34) 



71=0 



where 



CT,n+l(^) 

{ct,„+i(0) 



(35) 



cos(l];t)-z-Iisin(r!;t) 
+*<,n(0)?^ Sin in'j)}e^j>{-iA^t} 



A' 



T' 



+i4,„+i(0)^ sin {n'j)}exp{-tB^t} , 
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and 



A' 

a:. 

n 
Br, 



= 2(l-A) + (-2) = -2A, 

A , , A' 
= — (2n+l H , 



(36) 



= 2AV(^-?^)(^+l), 

^ A r„2 + 1)21 _ A! 



The purity of the central qubit (with an initial state 
(IT) + U))/V2) hi the symmetric phase (A > 1) can be 
determined as (see Appendix B) 



2 2 



2 1 2 

(^) (iV-n)(n+l)sin2(17;) . 



(37) 

In Case I (A' = A/iV), the purity P JSTl) can be further 
simplified as 



Especially, when the coupling strength between the cen- 
tral qubit and the bath is inversely proportional to the 
spin number of the bath, the central spin and the bath 
becomes decoupled, and the central qubit preserves its 
quantum coherence all the time. The present study not 
only demonstrates how the QPT influence the quantum 
coherence of the central qubit, but also establishes the 
connection between the LMG model and anti J-C model. 
In addition, our investigation may propose a new scenario 
to preserve quantum coherence of a central qubit in ex- 
perimental implementation of quantum computation. 
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APPENDIX A: DYNAMICS OF PURITY IN THE 
SYMMETRY BROKEN PHASE (0 < A < 1) 



P 



1 1 



1 



1 



A2 



sin^ {\t) 



(38) 



We see from Eq. ([55]) that P varies periodically 
and is independent of N though the coupling strength 
-2XWN - Sd in Eq. ([M]), as we have seen in Fig. 1. 
The physics behind Eq. ((38)) is that the mean "photon" 
number of the ground state is also A^-dependent, which 
countervails with the A'^-dependent coupling strength, 
leading to the A^-independent dynamical behavior of the 
purity P in Case I. 

On the other hand, in Case II {\' — \/ VA), the purity 
P ((57|) can be further simplifled as 



P = 



i + i{l-sin2 [v/Af(A2-l)t]} 



(39) 



Hence, the behavior of purity would not reach a steady 
pattern when A^ increases, as observed in Figs. 3 and 
4. The A^-dependence of the purity P in Case II stens 
from the coupling strength — 2A'\/ A — d'^d and the A^- 
dependent mean "photon" number of the ground state; 
the A^-dependence of them cannot countervail with each 
other. 



VI. SUMMARY 

We have studied the dynamic property of a central 
qubit coupled to an isotropic Lipkin-Meshkov-Glick bath. 
Two different types of coupling strength between the cen- 
tral qubit and the bath are considered. In both cases, the 
QPT of the bath is well revealed by the dynamic behav- 
ior of the central qubit. We have found that our model 
is equivalent to an anti J-C model under H-P transfor- 
mation when the bath is in the symmetry broken phase. 



For the anti J-C Hamiltonian ([24]) . the time evolution 
of a initial state |*(0)) = (|t) + |i))/\/2(g) |0) can be 
expressed as 



ik /A 

sin(r2oi)}exp i—t ] |t) ® |1)(A1) 



V2r2, 



iA 



COS [VLqI) + — — sin (f^ot) 



X exp 



||)^|O) + lIl0|O). 



The reduced density matrix [t) (Eq. (fT3l) ) of the system 
is then found to be 



TrB|*(t))(*(t)| 
2 



(A2) 



1 
2 

1 

+ 2 
1 

+ 2 
+h.c 



sin^ (flot) 



A 



IT) (Tl 



cos" (riot) + sin" (Qat) 



2nc 



cos (flot) - sin {Qot) 

Zi £0 



u)ai 



A 



exp ( ] IT) (il 



Applying Eq. (HH), we obtain the purity P (Eq.lgHl)) of 
the central qubit 



P = - 
4 



H 2 



1 + cos" (not) + ( ^ ) sin" (not) 



(A3) 



sin^ (r!o<) + 2 sin2(f7oi) 
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APPENDIX B: DYNAMICS OF PURITY IN THE 
SYMMETRIC PHASE (A > 1) 



In the large N limit, the reduced density matrix p^{t) 
fEq. lfTSll ) of the central qubit is derived as 



For a generalized anti J-C Hamiltonian (Eq. ([33ll ). the 
time evolution of an initial state |$(0)) = (|T) + U))/V2(g) 
|n) can be expressed as 



where 



V2{ 



i^?-^sin(r2;,_i<) 



(Bl) 



exp(-iAi_it) It) ® \n) + 



r^sin(r!:,_ii) 



exp(-iB„_ii) |i> ® \n-l) 

exp{~iB„t) \i) (g) \n) H 
exp{-iA„t) IT) «) |n+ 1) , 



jcos {n'^_,t) + sin {n:^^,t) 



TrB|<i>W) {m\ 

^IT) (Tl + ^li) ai + iexpH(A„-S„)t] X 



(B2) 



cos (r2;_it) - — 

x|T) {l\ + h.c. 



A(2n+1) + — 



Using Eq. ([Tl]), we obtain the purity P (Eq.([37l)) of 
the central qubit 



2A' 



(B3) 
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